Search for Invisible Decays of the J/V' in i/}{2S) — > tt+tt J/V' 
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Using iIj{2S) — > -r^-r' J/ijj events in a sample of 14.0 million ip{2S) decays collected with the 
BES-II detector, a search for the decay of the J/ip to invisible final states is performed. The J/ip 
peak in the distribution of masses recoiling against the tt+tt" is used to tag J/ip invisible decays. 
No signal is found, and an upper limit at the 90% confidence level is determined to be 1.2 x lO"'^ 
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for the ratio k(j/2S^+ -\ ■ This is the first search for J/xp decays to invisible final states. 



PACS numbers: 13.25.Gv, 95.30.Cq 

Invisible decays of quarkonium states such as J /i}) 
and T, etc., offer a window into what may lie beyond 
the standard model (SM) P, 0, [1]. This is because, 
aside from neutrinos, the SM includes no other invis- 
ible particles that these states can decay into. In the 
SM, the predicted branching fraction for J /ijj — )■ i/U is 

S(J/V' w) = 4.54 X 10"^ X B{Jlij ^ e+e") (1) 

with a small uncertainty (2-3%) [J, 3- However, new 
physics beyond the SM might enhance the branching 
fraction of J /ijj invisible decays. One possibility is the 
decay into light dark matter (LDM) particles medi- 
ated by a new, electrically neutral spin-1 gauge boson 
U ^ which could significantly increase the invisible de- 
cay rate [5]. On the other hand, C-even operators 
qq or 9750 do not contribute to invisible decays of 
the J/V- UM- 

Thus, there would be no contributions 
from a possible scalar or pseudoscalar exchange to the 
decay of J/?/" (3|jl6j]. 

Astronomical observations of a bright 511 keV 7- 
ray line from the galactic bulge have been reported 
by the SPI spectrometer on the International Gamma- 
Ray Astrophysics Lab (INTEGRAL) satellite 0. The 
corresponding galactic positron flux, as well as the 
smooth symmetric morphology of the 511 keV emis- 
sion, may be interpreted as originating from the anni- 
hilation of LDM particles into e+e" pairs It is of 
interest to search for such light invisible particles in 
collider experiments. CLEO reported an upper bound 
on T(IS') — > 7 -|- invisible, which is sensitive to dark 
matter candidates fighter than about 3 GeV/c^ and 
also provides an upper limit on the axial coupling of a 
new U boson to the h quark [8] . Recently, both CLEO 
and Belle reported upper limits on T(IS') invisible 
decays Q. The first experimental limits on invisible 
decays of the 77 and 77' mesons have been reported by 
the BES Collaboration [l^] ; these can be used to con- 
strain the mass of the LDM particles and couplings of 
new bosons to light quarks. Here, we present the first 
search for J/i/; decays to invisible final states. 

BES-II is the upgraded version of the BES-I detec- 
tor A 12-layer vertex chamber (VC) surrounding 
the beam pipe provides trigger and position informa- 
tion. This detector efficiently detects the presence of 
charged tracks over 97% of the total solid angle. A 
forty-layer main drift chamber (MDC), located radi- 
ally outside the VC, provides trajectory and energy 
loss {dE/dx) information for charged tracks over 85% 
of the total solid angle. The momentum resolution is 
CTp/p = 0.017 VTTp2 {p in GeV/c), and the dE/dx 
resolution for hadron tracks is ~ 8%. An array of 48 
scintillation counters surrounding the MDC measures 
the time-of-flight (TOE) of charged tracks with a res- 
olution of ~ 200 ps for hadrons. Outside of the TOE 
counters is a 12-radiation-length barrel shower counter 



(BSC) composed of gas tubes interleaved with lead 
sheets. This measures the energies of electrons and 
photons over ~ 80% of the total solid angle with an en- 
ergy resolution oiaE/E = 22%/ Ve {E in GeV). Out- 
side of the solenoidal coil, which provides a 0.4 Tesla 
magnetic field over the tracking volume, is an iron flux 
return that is instrumented with three double layers 
of counters that identify muons of momenta greater 
than 0.5 GeV/c. 

In order to detect invisible J/ip decays, we use 
ip{2S) TT^TT~J/ip and infer the presence of the J/V' 
from the J/ip peak in the distribution of mass recoil- 
ing against the tt+tt^. In this analysis, we use a 19.72 
pb~^ data sample collected at the peak of the tp{2S) 
resonance and a 6.42 pb^^ data sample collected off- 
resonance at a center-of-mass energy of ^/i = 3.65 
GeV. The data were recorded in the BES-II detec- 
tor operating at the Beijing Electron Positron Collider 
(BEPC). The data sample contains 14.0 million ij{2S) 
decays. 

In the search for J/ip invisible decays, invisible 
means nothing besides the two pions is seen in either 
the tracking or calorimetry systems of the detector. 
The charged-track trigger in BES-II requires at least 
one hit in the 48 barrel TOE counter array, one track 
in the VC and MDC^nd at least 100 MeV of energy 
deposit in the BSC [ll|, [l3] . This trigger is sensitive 
to two soft pions from the iIj{2S) decay of the signal 
events. 

Events with only two charged tracks with zero net 
charge are selected. Each charged track is required to 
be well fitted by a helix and to have a polar angle, 9, 
within the fiducial region | cos 6*1 < 0.8. To ensure that 
the tracks originate from the interaction region, we 

require Vxy = ^Jvf^-V^ < 2 cm and \Vz\ < 20 cm, 

where Vx,Vy, and Vz are the x,y, and z coordinates 
of the point of closest approach of each charged track 
to the beam axis. In addition, the momentum of each 
charged track must be less than 0.4 GeV/c. Particle 
identification (PID) is performed with the combined 
TOE and dE/dx information, and both charged tracks 
must be identified as pions. The invariant mass of the 
7r+7r~ pair is required to be larger than 0.35 GeV/c^. 

In BES-II, high momentum charged tracks with 
I cos 6*1 < 0.8 traverse the full radial extent of the MDC 
are detected and reconstructed efficiently. Tracks with 
lower momentum and/or at larger angles miss some 
or all of the MDC layers and may not be reliably re- 
constructed. However, tracks that originate from the 
center of the interaction point with a polar angle in 
the range | cos 9\ < 0.97 penetrate at least six layers of 
the inner vertex chamber; Although these tracks are 
not fully reconstructed, their presence can be inferred 
with high efficiency from track segments reconstructed 
from three or more VC hits. Events with indications 
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FIG. 1: Distribution of masses recoiling against the tt^tt" 
for V'(25) ■k'^tt~ J/^, J/tp invisible candidate events. 
Dots with error bars denote data. The lower histogram is 
the signal shape from tp{2S) tt^tt' J/if), J/tp — >■ f.i'^ j-i^ 
events. The upper histogram and the dashed curve are 
from the fit. 



of additional charged tracks anywhere in the region of 
angular coverage of the vertex chamber are rejected. 

Each reconstructed BSC cluster is required to have 
an energy greater than 10 MeV and to have a cluster 
profile consistent with that of a shower in the BSC. 
The number of unassociated neutral clusters, which 
do not match with either charged track in the MDC, 
is required to be zero in order to suppress backgrounds 
from J/ip decaying into neutral final states. 

Figure [T] shows the distribution of masses recoil- 
ing against the tt+tt" pair for candidate events. A 
clear J/tp peak is evident for the data taken at the 
'0(25'), while the smooth background under the J/tp 
peak from the QED contribution is consistent with 
the distribution obtained with off-resonance data at 
i/i = 3.65 GeV after normalizing to the luminosity at 

By measuring the ratio of invisible J/tp decays to 
J/tp — !• z^+zi", many uncertainties cancel, including 
those related to the number of J/tp decays, the soft 
pion tracking efficiencies, and the zero neutral clus- 
ter requirement. First, two soft pions are chosen 
with the same selection criteria as used for tp{2S) 
Tr~^TT~ J/tp, J/tp —> invisible. Then, we require that 
there are two muons in addition to the two soft pi- 
ons. Two selection criteria are used to identify muons. 
One is that the momentum of each track is larger than 
0.7 GeV/c. The other is that R, which is defined as 

be larger than 1.0 in order to reject J/tp e+e" 
events. Here, Egci and Esc2 denote the deposited en- 
ergies in the BSC, and pi and P2 denote their mo- 
menta. Figure [5] (a) shows the R distributions for 
data and Monte Carlo (MC) simulated events. For 



each event, the total energy of the four charged tracks, 
Etot, is required to be larger than 3.6 GeV, as shown 
in Fig. [2] (b). The small difference between the Etot 
distributions for data and MC simulation is caused 
by imperfect simulation of soft pions, an effect that 
is considered in the systematic error study. The re- 
quirement on the number of unassociated neutral clus- 
ters or extra charged tracks is the same as that for 
tp{2S) 7r+7r~ J/-0, J/tp invisible so that the sys- 
tematic uncertainties cancel in the ratio. Figure [3] 
shows the M^f^^3 distribution for 7r"'"7r~/x"'"/z~ final 
states for data and MC simulation after the above se- 
lection. Compared to the data distribution, the MC 
sample has a shift of 1.0 MeV/c^, which is also caused 
by the imperfect simulation of soft pions. The number 
of tp{2S) — > n'^TT^ J/tp, J/tp events selected 

is 43429 ± 208 where the error is statistical, and the 
corresponding MC-determined detection efficiency is 
14.5%. 
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FIG. 2: (a) Distributions of R for the two non-pion tracks 
in tp{2S) TT+TT- J/tp, J/tp ~* £+r. (b) Distributions 
of Etot for tp{2S) TT+Tv- J/tp, J/tp tJ.'^tJ.'- The his- 
togram denotes MC and dots with error bars denote data. 

Exclusive channels that are potential backgrounds 
are studied using full MC simulations in order to de- 
termine their contamination. The sources of the back- 
grounds can be divided into two classes. Class I is 




W^t+^Jrecoiling) (GeV/c^) 



FIG. 3: Distribution of mass recoiling against the tt 7r~ 
pair for tp{2S) tt+tt^p^/x" candidate events. The his- 
togram denotes MC simulation and dots with errors denote 
data. 
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mainly from final states without a J/ip^ such as the 
continuum processes e^e^ — > qq and r+r^ [q = u, 
d and s quarks). The smooth backgrounds under the 
J/ip peak in Fig. [T] are from "Class I" . Class II is from 
^j{2S) TT^TT^J/ip, where the J/tp decays into other 
modes than the true invisible final state. For this case, 
the J/ '0 decay products are either outside of the detec- 
tor acceptance, or are inside but are undetected. For 
these kinds of backgrounds, there is a, J/ip peak in the 
jyjrecoil distribution that is identical to the invisible 
signal. The largest sources of "peaking backgrounds" 
are from the J/ip decays to £^£~{£ = e,ii), where the 
two tracks tend to be back to back, so when one track 
escapes into the forward insensitive region, the other 
track tends to escape into the backward insensitive 
region. This is a geometric effect and its size is deter- 
mined from Monte Carlo simulation. Other sources of 
backgrounds are from decays to nn, where the neu- 
tron and anti-neutron leave no hit information in the 
detector. An additional potential background comes 
from J/ip decays to nhir'^, in which neither of the two 
photons from 7r° decay are within the detector accep- 
tance and the neutron and antineutron either have no 
hit information in the detector or miss the sensitive 
region of the detector. 

The systematic uncertainties on the number of ex- 
pected background events caused by the selection cri- 
teria are estimated using special event samples. For 
the dominant backgrounds from J/ tjj decays to i'^i^ , 
the uncertainty is mainly from the difference between 
data and MC simulation of the contamination from 
two missing leptons. We determine this systematic 
error using J/ip — s- KsK^tt^ , in which only the 
Ks and the charged pion are required in event selec- 
tion. This provides a large sample of (~38.6K) tagged 
charged kaons in the detector. With these, we mea- 
sure the probability for a charged kaon to be "invis- 
ible," i.e., to evade the extra track requirements, to 
be (3.49 ± 0.09)%. A MC simulation of this process 
gives an invisibility probability of (3.13±0.07)%. The 
difference of ~ 10% between the measured and sim- 
ulated probabilities is taken as a systematic error for 
peaking backgrounds from J/ip decays to and 
pp. 

In estimating the expected background contamina- 
tion from J/tp —f nn, we assume the branching frac- 
tion and its corresponding uncertainty for J/^p — s- nh 
are the same as those for J/ip — s- pp as is expected 
according to SU(2) symmetry. Similarly, since the 
measurements of branching fraction for J/tp nnir'^ 
is not currently available, we assume the branching 
fraction and its corresponding uncertainty for J/ip — > 
nnTr" are the same as those for J/ip — ppir^ . For the 
J jip ^ nn background, the main uncertainty is from 
the difference between data and MC simulation in the 
nn rejection rate based on the requirement of zero 
clusters in the BSC. The distribution of the number of 
clusters associated with nn is obtained from two con- 
trol samples, J/?/; — > pfvK~ +c.c. and J/ip —^ ■n'^-K~pp. 
First, the distribution of the number of clusters ex- 



TABLE I: Expected number of events (Ntg) and efficien- 
cies for peaking backgrounds. 



Background channel 


Efficiency 


Expected 


{lp{2S) TT+TT-J/i:, J/ij 


(%) 






0.964 


2543 ± 254 




0.907 


2393 ± 240 


nn 


10.46 


1011 ±85 


PP 


0.434 


42 ± 13 


nnn'^ 


0.486 


29± 10 


Total 




6018 ±514 



pected for ppTT^n^nn can be obtained by combining 
the pnTr~ and pmr^ samples. Second, using the con- 
trol sample, J/ip 'k~^tt~pP, the cluster distribution 
for nn can be obtained by subtracting the pp'K^-K~ 
distribution from the ppTr~^T:~nn distribution. The 
distributions of the number of clusters between data 
from the control sample and MC simulation are com- 
pared, and the difference between MC and data for the 
number of nn events with zero BSC clusters is a 7.6% 
effect, which is taken as a systematic error. The es- 
timated contributions from the peaking backgrounds 
are summarized in Table HI The error on the expected 
number of events for each peaking background in- 
cludes the contributions from both the branching frac- 
tion uncertainty and the estimated systematic error. 
For J/ip £^i^ and pp, the dominant uncertainty 
arises from the geometrical acceptance and their as- 
sociated errors are entirely correlated. 

A fit is used to extract the number of J/ip 
events in the distribution of mass recoiling against 
the TT+n- in the range 3.0 GeV/c^ < M^^^oil ^ 3 2 
GeV/c^. Here, the shape used to describe the sig- 
nal comes from the it^it~ recoil mass spectrum from 
the i^{2S) TT^TT~J/tp, J/tp ^ control sample, 

and not from the MC simulation due to differences 
between the data and MC simulation (see Fig. [3]). 
The Class I background is represented by a second- 
order polynomial. The fit, with /ndf = 62.6/75 
and shown in Fig. [H yields 6424 ± 137 events in the 
peak, which includes the contributions from both sig- 
nal and peaking backgrounds, since they have the 
same probability density functions (PDF) in the fit. 
After subtracting the expected backgrounds listed in 
Table U from the fitted yields, we get the number of 
ip{2S) Ti^TT' J /"ip, J /ip invisible events to be 
406 ± 532. 

The estimated uncertainties that do not cancel in 
the ratio are described here. The systematic uncer- 
tainty caused by the tracking efficiency in the MDC 
for the two muons in the control sample is estimated 
to be 4% [1^. The systematic error for the require- 
ment on Etot is 1.7%, and the uncertainty caused by 
the R requirement in the selection of the control sam- 
ple, ip{2S) TT+TT^ J/i/j, J/tp ^ A'^M I is determined 
to be 1.0% by a comparison with and without the R 
requirement. The uncertainty from the background 
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shape, which is used in the fit to M^^'^'^}^ for Class 
I, is found to be neghgible. The uncertainty associ- 
ated with bin-size and the range of the fit is 1.3%. 
In order to investigate the uncertainty caused by the 
trigger, we check the four trigger channels used in 
the BES-II experiment. The largest systematic er- 
ror comes from uncertainties in modeling the energy 
threshold requirement in the BSC. A conservative es- 
timate, based on MC simulation, is that this corre- 
sponds to at most a 1% uncertainty in the detection ef- 
ficiency. The uncertainty of B{J/'ip — *■ is taken 
from the PDG [U- The total uncertainty, which is 
determined by the sum of all sources in quadrature, 
is 4.9%. Taking this systematic uncertainty into ac- 
count, the upper limit for the number of events of 
■0(25') TT+vr-J/i/', J/'0 invisible is 1285 at the 
90% confidence level, or a central value of 406^^333 at 
the 68.3% confidence level from the Feldman-Cousins 
frequentist approach [l5l |. 

Finally, the upper limit on the ratio of B{J/'iJj — > 
invisible) to B{J/ip — > is determined from the 

relation 



B{J/ip — > invisible) 



N, 



< 



J/V 

UL 



/^invisible 



1.2 X 10~ 



(2) 



where iV^{f (1285) is the 90% confidence-level upper 
limit on the number of iJj{2S) — > n'^ Tr~ J / tp , J / — s- 
invisible events, £invisible('^^-^'^) MC- 
determined signal efficiency, (43429 ± 208) is 

the mmibcr of events for ip(2S) — > tt^ J / ijj , J/ijj — > 
and e^i'^- (14.5%) is the detection efficiency for 
that decay mode. In addition, the two-sided interval 
of the number of the measured events for invisible de- 
cays, iV;^'ll»6ie' is (73, 945) at the 68.3% confidence 



level. The corresponding two-sided interval of the ra- 
tio "'^itii^fif^' is (0.66 X 10-^ 8.6 X 10-^). 

In summary, we performed the first search for invisi- 
ble decays of the J/ip using tp{2S) ■n^ir^ J/ip events 
detected in a sample of 14.0 million tp{2S) decays. The 

upper limit on the ratio ^^^(^//^"^J+^^l'^^ at the 90% 
confidence level is 1.2 x 10~^. This measurement im- 
proves by a factor of 3.5 the bound on the product of 
the coupling of the U boson to the c quark and LDM 
particles as described in Eqs. (25) and (26) of Ref. 0]. 
One now has, for a Majorana LDM particle x as in 
Eq. (26) of Ref. 0, a limit of \c^fcv\ < 8.5 x IQ-^, 
which is almost a factor of 2 stronger than the cor- 
responding limit |cpj-/f,v| < 1.4 x 10~^ derived from 
the invisible decays of the T(15') as described in Eq. 
(106) in Ref. Q, where and fcv ihv) denote the 
U boson couplings to the LDM particle x and c (&) 
quark. We expect a more precise measurement can be 
obtained in the future BES-III experiment. 
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